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ABSTRACT 

We study two dimensional turbulent magnetic reconnection in a compressible fluid in the gas pressure dominated limit. We use 
open boundary conditions and start from a Harris current sheet configuration with a uniform total pressure. A small perturbation to 
the vector potential initiates laminar reconnection at the Sweet-Parker rate, which is allowed to evolve for several dynamical times. 
Subsequently sub-Alfvenic turbulence is produced through random forcing at small wave numbers. The magnetic field topology 
near the current sheet is strongly affected by the turbulence. However, we find that the resulting reconnection speed depends on the 
resistivity. In contrast to previous results in three dimensions, we find no evidence for fast reconnection. The reconnection speed 
exhibits large variations but the time averages increase smoothly with the strength of the turbulence. 

Key words, galaxies: magnetic fields — physical processes: MHD — physical processes: turbulence — methods: numerical 



1. Introduction 

Radio observations of a wide variety of astrophysical bodies (the 
Sun, spiral galaxies, the Earth) show that a magnetic field is 
usually present (see [Priest & F orbes 2000 for review). These 
magnetic fields have significant large scale components, i.e on 
the scale of the objects themselves. The origin of these fields is 
typically ascribed to the operation of a large scale magnetohy- 
drodynamic dynamo (Ruzmaikin et al. '1988': 'Parker 1979 1992 



[Hanasz et al.|,2 004; Lazarian & Vishniac 20081. Since dynamo 
theory involves the twisting and folding of field lines it is im- 
portant that there is some process which can lead to efficient 
smoothing of the small scale components of the field. In other 
words, we need to invoke some kind of fast local magnetic dif- 
fusion. In spite of the fact that astrophysical fluids are turbulent 
(see Armstrong et al. 1995 ; Horbury & Balogh 2001 ; Elmegreen] 
^Scalo„2004, ,McKee & Ostriker,2007j), the concept of turbu- 
lent magnetic difFusivity (see Blackman & Field|2008) l, a popu- 



ill founded ( 


Cattaneo & Vainshtein 1991 , Gruzinov & Diamond 


|1994 Vains 


itein & Cattaneo 


1992|i. In particular, it does not ad- 



change their topology. In an ionized plasma, Ohmic difFusivity 
fails by many orders of magnitude to supply the required mag- 
netic diffusion. 

In order for astrophysical dynamos to function smoothly 
there must be a process which allows reconnection to proceed 
at speeds characteristic of local dynamical velocities. Since rms 
fluid velocities are often comparable to the local Alfven speed, 
in practice this requirement is indistinguishable from Vrec ~ Va- 
This is called fast reconnection, meaning that it does not depend 
on resistivity or depends on the resistivity logaiithmically (see 



from studies of solar flares ( 


Yokoyama & Shibata|1995 Innes et 


al.|1997 Pagano et al.|2008 





The idea that there could be some way of producing fast 
magnetic reconnection even in highly conducting fluids is not 
new (Moffat 1978;I Krause & Ra dler 1980), but early models of 
reconnection (see Parker"1957; Sweet 1958, e.g.) using realistic 
astrophysical temperatures and densities gave a very slow recon- 
nection rate, ~ VaS^^^^, where S = jj/CVaL) is the Lundquist 
number, rj is the resistivity, and L is the size of the cuiTent 
sheet. It was |Petsch ek ( 1964) who for the first time introduced a 
model for fast reconnection with a rate proportional to (log 5) ' . 
Subsequent numerical simulations and theoretical analyses have 
shown that the Petschek reconnection rate is only attainable in 
very restricted circumstances. For instance, a modified version 
can stably persist in a collisionless plasma (see [Drake et al.| 
2006a e.g.). This means that the length of the current sheet 
should not exceed approximately 50 electron mean free paths 
(Uzdensky 2006 Yamada et al. 2006[ l. This condition cannot 
be satisfied in many astrophysical environments, e.g. in the in- 
terstellar medium (Vishni ac & Lazarian 1999). In a collisional 
plasma the X-point region required for Petschek reconnection 
will collapse to the Sweet-Parker geometry for large S (Biskamp 
[T996[ ). 

The failure of the Petschek model has increased interest in 
the role of turbulence in reconnection. This interest has been 
further stimulated by the fact that the turbulence is ubiquitous in 
astrophysical environments where reconnection occurs, e.g. the 
ISM, stars, the Sun and accretion disks ( [Ruzmaikin et al.|1988[ l. 
The idea that turbulence can affect reconnection has a long his- 
tory, although usually studied in two-dimensions (2D) ( ,Priest & 
Forbes |2000[ l. Several researchers have approached this prob 



lem numerically, e.g.[Matthaeus & Lamkin'(1985"1986');'Fan et| 
al. (2004 200^; ISemdio et al.[ (2009); Loureiro et al. (200917 
They found that it was possible to get many features expected 
from reconnection theory, i.e. large and small-scale magnetic is- 
lands, fluid jetting, current filamentation and that the maximum 
reconnection speed was higher for more powerful turbulence and 
exceeded the Sweet-Parker rate (Fan et al.f2004 2005). 
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The most interesting result of their calculations (Fan et al. 
[2004[[2005j ) was that under the influence of turbulence they ob- 
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served a two-step process of magnetic reconnection: beginning 
with a slow Sweet-Parker mechanism and changing later to a 
faster reconnection state that they identified with the Petscheck 
process. In a similar manner, the most recent paper ( [Loureiro] 
|et al.|[2009) also observed that the presence of turbulence sig- 
nificantly enhances the reconnection rate. They also found that 
for a given value of diffusion and above a critical value of the 
turbulent injection power the reconnection process accelerates 
substantially. However, in all these simulations the authors used 
periodic boundary conditions, which prevents inflow or outflow 
and conserves the total magnetic flux. This made the actual re- 
connection rates difficult to evaluate. In particular, it was impos- 
sible to calculate the average reconnection rate, which is impor- 
tant in view of the large fluctuations induced by turbulence. In 
addition, one can argue that the reconnection rate in these sim- 
ulations is influenced by the conservation of the total magnetic 
helicity (Blackman'2000). These difficulties can be avoided by 
using boundary conditions that allow the inflow and outflow of 
plasma and magnetic flux ( Kowal et al.|2009l l. 

More recently interest in the magnetic reconnection process 
has moved from 2D magnetic configurations to more realistic 
and generic three-dimensional ones. In particular, |Lazarian"&] 
VishniacI ([T999| (LV99) and [Lazarian & Vishniac|pOOO I) pro- 



posed that in three dimensions a stochastic magnetic field com- 
ponent can dramatically enhance reconnection rates, leading to 
reconnection speeds comparable to the local turbulent velocity. 
Their model is based on the Sweet-Parker reconnection scheme, 
with a long narrow current sheet between two regions of dra- 
matically diff'erent polarizations but similarly strong magnetic 
fields, but includes the effects of turbulence and substructure in 
the magnetic field. This has two principal effects. First, in three 
dimensions many independent patches of magnetic field come 
into contact with the current sheet and undergo reconnection. 
Second, the outflow of plasma and shared magnetic flux happens 
not over a microscopically narrow region determined by Ohmic 
diffusion, but through a substantially wider region determined 
by field wandering. Neither effect is present in two dimensions, 
although the formation of magnetic islands in two dimensions is 
roughly similar to the broadening of the outflow. Together these 
effects are sufficient to trigger fast magnetic reconnection. In this 
model the reconnection rate does not depend on the Ohmic re- 
sistivity, but is determined only by turbulence, in particular by 
its strength and injection scale. This fast reconnection model 
has been tested numerically by Kowal et al. 12009) using in- 



flow/outflow boundary conditions and using a wide range of in- 
jection scales and power for the turbulence. These simulations 
have confirmed all of the predicted features of the LV99 model, 
including the insensitivity to the Lundquist number. 

In this paper we return to two dimensional reconnection us- 
ing the same inflow/outflow boundary conditions. We have two 
objectives in this work. First, since the explanation for fast re- 
connection advanced in LV99 was intrinsically three dimen- 
sional, we are interested in examining the effects of dimension- 
ality on the reconnection rate. Second, this work constitutes an 
examination of the importance of boundary conditions on the 
two dimensional model and a test of claims for fast reconnection 
in two dimensions. Our numerical model of turbulent reconnec- 
tion in the ISM is calculated in a 2D box with open boundary 
conditions. We use a Harris current sheet setup as an initial con- 
figuration. Reconnection develops as a result of an initial vector 
potential perturbation. We do not drive reconnection. We solve 
2D non-ideal normalized isothermal MHD equations numeri- 
cally while varying the power and scale of turbulence, and the 
magnetic resistivity ( [Kowal et al.|2009| . 



The plan of this paper is as follows. In Sect. |2]we describe 
the numerical setup and input parameters, and we present our 
method to measure the reconnection rate. The results are dis- 
cussed in Sect. [3] where we analyse the time evolution of mod- 
els with different power of injecting turbulence. We present the 
dependencies of the reconnection rate on: power of turbulence, 
injection scale, viscosity as well as uniform and anomalous re- 
sistivities. We also check the influence of initial magnetic field 
configuration, boundary conditions and method of driving turbu- 
lence. We discuss our results in Sect.|4]and give our conclusions 
in Sect.|5] 



2. Methods 

2.1. Basic Equations 

We study the problem of magnetic reconnection in the presence 
of weak turbulence using the magnetized fluid approximation 
governed by the isothermal non-ideal MHD equations of the 
form: 

dp 



dt 



H- V ■ (pv) = 0, 



dpv 

Ik 



+ V ■ (pvv + PJ - BB) =/+ vAv, 



dA 

~dt 



V xB - rjV xB, 



(1) 



(2) 



(3) 



where p is the gas density, v is the fluid velocity, A is the 
magnetic vector potential, B = V x A is the magnetic field, 
Pt = c]p+B^ l^n is the total pressure, / is the Kronecker delta, Cj 
is the isothermal speed of sound, t] is the resistivity coefficient, 
V is the viscosity, f - pa represents the forcing term, and a is a 
random acceleration. 

We solve the MHD equations using the same code as in 
Kowal et al.|d2009) based on the following methods: a higher- 



order shock-capturing Godunov-type scheme, the essentially 



non oscillatory (ENO) spatial reconstruction (see Londrillo & 
[Del Zanna|[20D0i [Del Zanna eraL] [2003 ), a multi-state Harten- 
Lax-van Leer (HLLD) approximate Riemann solver for isother- 
mal MHD equations ('Mignone'2007 ) and Runge-Kutta (RK) 
time integration (see )Del Zanna et al.^2003 , e.g.). The choice of 
HLLD Riemann solver guarantees a good solution for the Alfven 
wave propagation, which is important in our model, since most 
of the kinetic energy is transported through Alfven waves (see 
Kowal et al. 2009). The divergence of the magnetic field must 
be kept zero everywhere at all times (V ■ B = 0). To satisfy this 
condition the field interpolated constraint transport (CT) scheme 



based on the staggered grid is used (see Londrillo & Del Zanna 
2000, e.g.). 

Some selected simulations that we perform include anoma- 
lous resistivity modeled as 



Jcrit 



I H 



Jcrit 



(4) 



where rj,, and t/q describe uniform and anomalous resistivity co- 
efficients, respectively, jc-ii is the critical level of the absolute 
value of current density above which the anomalous effects start 
to work, and // is a step function. However, for most of our sim- 
ulations T]a = 0. 
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Fig. 1. 2D magnetic field configuration in our problem. The grey 
area describes the diffusion region where the incoming field lines 
reconnect. The longitudinal and transverse scales of the diffusion 
region are described by the parameters A and 6, respectively. We 
use inflow and outflow boundary conditions at X and Y direc- 
tions, respectively. 



2.2. Initial and Boundary Conditions 

We numerically investigate the turbulent magnetic reconnection 
model in 2D in a computational box of size x L,,, where 
= 1 and Ly - 2 with a spatial resolution of 1024x 2048 
grid zones in the x and y directions, respectively. Figure [T] 
shows a 2D visualization of the reconnection problem setup. 
The domain contains two regions of oppositely directed mag- 
netic lines separated by a diffusion region with a thickness 6 and 
a length A, where the magnetic lines reconnect and the prod- 
uct of this process is ejected along the X direction with a speed 
Voui (see Figure [T]). The initial magnetic field configuration is 
described by a Harris current sheet B^{x,y) - tanh(3'/0), 
where Bjo is the initial strength of the anti-parallel magnetic 
field component. The Sweet-Parker reconnection is triggered by 
a small initial perturbation of the vector potential 6A{x,y) - 
6Bocos{27rx)exp[-(y/d)^], where d and dBo denote the thick- 
ness of the perturbed region and the strength of the initial per- 
turbation, respectively. The initial setup is completed by setting 
the density profile from the condition of the uniform total pres- 
sure Pt{t = 0, x,y) = const and setting the initial velocity to zero 
everywhere. 

The code ( |Kowal et al.|2009| ) uses dimensionless units. The 
velocity and magnetic field are expressed in units of the char- 
acteristic Alfven speed Va = IB.iol/ -\/pq, where the initial anti- 
parallel component of magnetic field = 1 and the density 
Pq - 1 far from the diffusion region so that V/i = 1 far from 
the diffusion region in all models. Time is measured in units of 
the Alfven transit time = 1/Va- The speed of sound is set to 
Cj = 4 (the plasma /3 = p/ p,„„g is 32.0 for all models). We vary 
the resistivity coefficient ?;„ between values 3 ■ 10""^ and 5 ■ lO"-' 
(in dimensionless units). In the models which include anomalous 
effects, we vary the anomalous resistivity coefficient rja between 
0.0 and 3 ■ 10"^. The parameters describing the initial perturba- 
tion are set to SBq = 0.05 and li = 0. 1 . 

We use outflow boundary conditions along the X direction 
and inflow boundary conditions along the Y direction setting the 
normal derivatives of the fluid variables (density and momen- 
tum) to zero. In the treatment of vector potential A at the bound- 
ary we set its components transverse to the considered boundary 
using the first order extrapolation, while the normal derivative 



of the normal component is set to zero. This guarantees that all 
waves generated in the system are free to leave the box without 
significant reflections. For a more detailed description of these 
boundary conditions, including their advantages and drawbacks, 
we refer the reader to |Kowal et al.| ( |2009l l. 

2.3. Model of Turbulence 



In order to drive turbulence in our model we follow iKowal et aO 
( 2009| l and use the method proposed by 'Alvelius ( 1999)1. The 
forcing is driven across a specified distribution of wave vectors. 
Here we use a Gaussian profile around a shell in Fourier space 
with a radius which determines the injection scale /,„j. The forc- 
ing is random in time, and therefore is uncorrected with any of 
the time scales of the turbulent flow. For the same reason the 
power input is defined purely by the force-velocity correlation. 
The driving is completely solenoidal and does not directly pro- 
duce density fluctuations. 

Turbulence is introduced at a given injection scale and grows 
gradually in time until it reaches the desired amplitude corre- 
sponding to the turbulent power P,„j. We drive turbulence to 
its saturation level over one Alfvenic time from ti,eg = 9 to 
tend - 10. According to the LV99 model, the injection scale and 
turbulent power determine the rate of reconnection in 3D. Thus, 
in our model we test this correlation in 2D by changing these 
properties of turbulence. 



2.4. Input Parameters 

Our simulations can be divided into five groups. In each of them 
we analyze the dependence of the reconnection rate by changing 
only one of the crucial parameters: for models PD - the power of 
turbulence P,„j, for models SD - the injection wavenumber kf, 
for models VD - viscosity v, for models RD and AD the uniform 
T]u and anomalous rja resistivities, respectively. In Table[T]we list 
parameters of all the models presented in this paper. 

Among all parameters of the model we list those which vary, 
i.e. the uniform and anomalous resistivities, 77„ and rja, respec- 
tively, the uniform viscosity v, the power of turbulence Pi„j and 
its injection wavenumber ki„j with the half- thickness of the in- 
jection shell Akinj, the number of perturbed Fourier components 
of velocity Nf and the amplitude of perturbation v/ at the injec- 
tion scale. In addition, we include the mean velocity amplitude 
obtained at the end of each simulation. 



2.5. Reconnection Rate Measure 

In the case of laminar reconnection (Sweet-Parker and Petschek) 
its rate can be measured by averaging the inflow velocity y,„ 
divided by the Alfven speed Va over the inflow boundaries, i.e. 



<V/„/Va> 







Vy 






>-y,,m, 


Va 


J 



(5) 



Since we have two X boundaries, located al y - y^in and 
y - ymax^ we need to take half of the resulting integral. This 
measure works well for laminar reconnection, when the system 
is perfectly stable and where the time derivative of the magnetic 
flux is zero. In the presence of turbulence, however, this time 
derivative can fluctuate or the turbulence in the center of the box 
could affect the flow of the plasma. In this way we would get a 
flow of magnetic flux without the presence of reconnection. In 
order to include all effects contributing to the change of magnetic 
flux, we use a new more general measure of the reconnection rate 
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Table 1. List of models 



Name 


ri. [10-'] 


[10-^ 


V [10-"] 


p 

^ in } 






Nf 


Vf 


<M> 


PD 


0.5 


0.0 


0.0 


0.005 


12 


0.5 


68 


0.00006 


0.03983 




0.5 


0.0 


0.0 


0.01 


12 


0.5 


68 


0.00009 


0.04912 




0.5 


0.0 


0.0 


0.05 


12 


0.5 


68 


0.00019 


0.08401 




0.5 


0.0 


0.0 


0.1 


12 


0.5 


68 


0.00028 


0.10645 




0.5 


0.0 


0.0 


0.5 


12 


0.5 


68 


0.00062 


0.21685 




0.5 


0.0 


0.0 


1.0 


12 


0.5 


68 


0.00088 


0.28657 


SD 


0.5 


0.0 


0.0 


0.1 


8 


0.5 


48 


0.00033 


0.11618 




0.5 


0.0 


0.0 


0.1 


12 


0.5 


68 


0.00028 


0.10645 




0.5 


0.0 


0.0 


0.1 


16 


0.5 


112 


0.00029 


0.08927 




0.5 


0.0 


0.0 


0.1 


20 


0.5 


112 


0.00023 


0.06626 




0.5 


0.0 


0.0 


0.1 


24 


1.0 


36 


0.00042 


0.05534 




0.5 


0.0 


0.0 


0.1 


28 


1.0 


32 


0.00058 


0.04633 




0.5 


0.0 


0.0 


0.1 


32 


1.0 


32 


0.00056 


0.03997 


RD 


0.3 


0.0 


0.0 


0.01 


12 


1.0 


68 


0.00009 


0.05734 




0.5 


0.0 


0.0 


0.01 


12 


1.0 


68 


0.00009 


0.04912 




0.6 


0.0 


0.0 


0.01 


12 


1.0 


68 


0.00009 


0.04471 




0.7 


0.0 


0.0 


0.01 


12 


1.0 


68 


0.00009 


0.04889 




0.8 


0.0 


0.0 


0.01 


12 


1.0 


68 


0.00009 


0.05405 




0.9 


0.0 


0.0 


0.01 


12 


1.0 


68 


0.00009 


0.04813 




1.0 


0.0 


0.0 


0.01 


12 


1.0 


68 


0.00009 


0.05481 




2.0 


0.0 


0.0 


0.01 


12 


1.0 


68 


0.00009 


0.05542 




3.0 


0.0 


0.0 


0.01 


12 


1.0 


68 


0.00009 


0.06214 




4.0 


0.0 


0.0 


0.01 


12 


1.0 


68 


0.00009 


0.06774 




5.0 


0.0 


0.0 


0.01 


12 


1.0 


68 


0.00009 


0.07407 




0.3 


0.0 


0.0 


0.1 


12 


1.0 


68 


0.00028 


0.11583 




0.5 


0.0 


0.0 


0.1 


12 


1.0 


68 


0.00028 


0.10649 




0.6 


0.0 


0.0 


0.1 


12 


1.0 


68 


0.00028 


0.10763 




0.7 


0.0 


0.0 


0.1 


12 


1.0 


68 


0.00028 


0.10122 




0.8 


0.0 


0.0 


0.1 


12 


1.0 


68 


0.00028 


0.10221 




0.9 


0.0 


0.0 


0.1 


12 


1.0 


68 


0.00028 


0.10036 




1.0 


0.0 


0.0 


0.1 


12 


1.0 


68 


0.00028 


0.09854 




2.0 


0.0 


0.0 


0.1 


12 


1.0 


68 


0.00028 


0.08844 




3.0 


0.0 


0.0 


0.1 


12 


1.0 


68 


0.00028 


0.08075 




4.0 


0.0 


0.0 


0.1 


12 


1.0 


68 


0.00028 


0.07236 




5.0 


0.0 


0.0 


0.1 


12 


1.0 


68 


0.00028 


0.06987 


AD 


0.5 


0.0 


0.0 


0.1 


12 


1.0 


68 


0.00028 


0.10645 




0.5 


0.5 


0.0 


0.1 


12 


1.0 


68 


0.00028 


0.10534 




0.5 


1.0 


0.0 


0.1 


12 


1.0 


68 


0.00028 


0.11071 




0.5 


2.0 


0.0 


0.1 


12 


1.0 


68 


0.00028 


0.11252 




0.5 


3.0 


0.0 


0.1 


12 


1.0 


68 


0.00028 


0.11465 


VD 


0.5 


0.0 


0.08 


0.1 


12 


1.0 


68 


0.00028 


0.10372 




0.5 


0.0 


0.09 


0.1 


12 


1.0 


68 


0.00028 


0.09623 




0.5 


0.0 


0.10 


0.1 


12 


1.0 


68 


0.00028 


0.09041 




0.5 


0.0 


0.20 


0.1 


12 


1.0 


68 


0.00028 


0.08472 




0.5 


0.0 


0.30 


0.1 


12 


1.0 


68 


0.00028 


0.08366 




0.5 


0.0 


0.50 


0.1 


12 


1.0 


68 


0.00028 


0.07317 




0.5 


0.0 


0.60 


0.1 


12 


1.0 


68 


0.00028 


0.06957 




0.5 


0.0 


0.80 


0.1 


12 


1.0 


68 


0.00028 


0.06868 




0.5 


0.0 


1.00 


0.1 


12 


1.0 


68 


0.00028 


0.06019 




0.5 


0.0 


2.00 


0.1 


12 


1.0 


68 


0.00028 


0.04371 




0.5 


0.0 


3.00 


0.1 


12 


1.0 


68 


0.00028 


0.03852 




0.5 


0.0 


4.00 


0.1 


12 


1.0 


68 


0.00028 


0.03394 




0.5 


0.0 


5.00 


0.1 


12 


1.0 


68 


0.00028 


0.02912 



as described in |Kowal et al.| ( |2009] l, which in 2D is given by the 
following simplified equation: 



1 



-5, 



(sign(B,)£,)|,, - (sign(B,)£,) 



dy\B, 



(6) 



where |B,,,oo| is the asymptotic absolute value of and £, is 
the Z component of the electromotive force. A more complete 



discussion of this new reconnection measure can be found in 
IKowal etal.| ( |2009l ). 

3. Results 

3.1. Laminar Reconnection 

During the first stage of our simulations, before we start driv- 
ing turbulence, the system evolves to reach the stationary state 
of Sweet-Parker reconnection. After that the influence of turbu- 
lence on the evolution of our system can be studied precisely. As 
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Fig. 2. Time variations of the total mass M (solid line), magnetic 
Emag (dashed line) and kinetic Eki„ (dotted line) energies during 
the Sweet-Parker reconnection stage with the uniform resistiv- 
ity rju - 5 ■ 10"'*. For clarity the kinetic energy Ei^i„ has been 
amplified by a factor of 100. 




Time 



Fig. 3. Time variation evolution of the reconnection rate during 
the Sweet-Parker stage with the uniform resistivity t/,, = 5 ■ lO"'* 



mentioned above we start our simulations with a configuration 
of the oppositely directed initial magnetic fields with a small 
magnetic perturbation. The perturbation initiates Sweet-Parker 
reconnection which reaches a stationary state in a few Alfen time 
units which lasts until we start injecting turbulence. 

We assume that we get the steady state of Sweet-Parker re- 
connection when the total mass, reconnection rate, kinetic and 
magnetic energies show very small time derivatives. Figure [2] 
shows the evolution of mass, kinetic and magnetic energy dur- 
ing the laminar reconnection stage. In the beginning all these 
quantities change slightly but then reach almost constant values. 
The reconnection rate (Figure|3]l initially increases until t ~ 1.3, 
when it starts to oscillate, finally stabilizing after a few time 
units. 

Figure |4] shows the topology of the velocity (left panel) and 
magnetic fields (middle panel), and the absolute value of the 
current density (right panel) just before we start injecting tur- 
bulence. The brighter shades correspond to large values of the 
displayed quantities. The velocity field and magnetic field are 
shown in the form of texture. The initial oppositely directed 
magnetic field lines are transported to the middle of the box. 
The Y-component of magnetic field emerges from the current 
sheet and is ejected near the midplane through the left and right 
boundaries. The reconnection process in the diffusion region and 
the ejection of By cause the strong outflow of gas which is clearly 



0.00 




Fig. 5. The evolution of the reconnection speed for a model 
with Pinj — 0.1, a mean wavenumber of the injected turbulence 
kf - 12 and uniform resistivity rju - 5 ■ 10""*. Marked time steps 
correspond to those presented in Figure |6] 

visible in the left panel of Figure |4] The system reaches a steady 
state despite this loss of mass due to the inflow of gas through 
the top and bottom of the computational box. 

The absolute value of current density is shown in the right 
panel of Figure]?] As expected the highest values of the current 
density appear in the midplane and determine the diffusion re- 
gion. The state of Sweet-Parker reconnection described above is 
stationary and sufficient to study the influence of turbulence on 
the reconnection process. 

3.2. Reconnection in the Presence of Turbulence 



In this section we present the influence of turbulence on the re- 
connection process. Turbulence is injected at time f = 9 and 
at a given injection scale li„j cc kj^ in the vicinity of the mid- 
plane. We gradually increase the strength of turbulence during 
one Alfven time, thus at time f = 10 the power reaches its input 
value defined by P„y . 

Figure |5] shows the reconnection rate obtained for a model 
with uniform resistivity /;„ = 5 ■ 10""*, a turbulent power Pj„j = 
0.1 and an injection wavenumber kf = 12. Adding turbulence to 
the system results in slight fluctuations of the reconnection rate 
until f = 1 1 .5. Next, the reconnection rate increases significantly 
until it reaches the maximum at time f ~ 13 and then drops. We 
can distinguish four such maxima, which are roughly separated 
by two Alfven times (t ~ 13, f ~ 15, f ~ 17, f ~ 19). In Figure]5] 
we mark three time steps corresponding to the same time stages 
shown in Figure ]6] where we plot the current density, magnetic 
and gas velocity fields. 

In Figure |6] (top and middle row) we see that the structure 
of the magnetic and velocity fields are considerably different 
than in the case of Sweet-Parker reconnection (see Figure ]4|. 
At fi = 12.15 (Figure]6] first column, top panel) we are injecting 
turbulence with the maximum power f ,„y = 0. 1 in a large vol- 
ume surrounding the the midplane. Thus, in this region the ve- 
locity field is strongly perturbed and mixed. Although the topol- 
ogy of the velocity field is very complicated, we can distinguish 
the main direction of the velocity fluctuations, which is parallel 
to the mean magnetic field. 

A smaller number of distinct features in the velocity field 
are pointed perpendicular to the mean magnetic field. Close to 
the midplane the magnetic field lines change their directions and 
are substantially reduced in strength. Thus, velocity fluctuations 
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Fig. 4. Absolute value of the current density distribution (right panel), the magnetic (middle panel) and velocity field (left panel) 
topology during the Sweet-Parker stage for time f = 9 in a model with the uniform resistivity rj,, = 5 ■ 10""^. Brighter shades 
correspond to the larger values of the displayed quantities. All panels show the moment just before we start injecting turbulence. 



can bend magnetic field lines in this region. Another noticeable 
difference in comparison to the Sweet-Parker configuration is 
a significant change of the current density distribution. Adding 
turbulence to the system creates a very complex configuration 
of the magnetic field, so that we observe multiple reconnection 
events happening at the same time (Figure|6] fi = 12.15, bottom 
panel). 

The magnetic field configuration looks quite different at the 
next time step t2 - 13.85 (Figure |6] second column, middle 
panel). In the vicinity of the midplane we observe the formation 
of a magnetic island. Also the velocity field (Figure |6] second 
column, top panel) is more mixed in this region. However, the 
outflow velocity of the gas velocity is quite low, which is con- 
firmed by the small reconnection rate (Figure |5]l. 

As the simulation proceeds a large loop of magnetic field 
lines moves to the left boundary (Figure |6j fa - 15.15, last col- 
umn). This causes a powerful outflow of gas, which results in 
a violent and rapid growth of the reconnection rate (Figure |5]l. 
We see that the maxima of the velocity field associated with this 
loop are where the value of magnetic field is low. The area of the 
magnetic island f2 - 13.85 andfs = 15.15 is also well defined by 
the current density distribution (Figure |6] bottom panel, middle 
and right column). 

Figure |7] shows the evolution of the reconnection speed for 
Pinj = 0.5, five times larger than in the previous case. As before 
the reconnection rate does not reach a stationary state. We ob- 
serve a few maxima, however they do not appear periodically as 
in the previous case (Figure |5]). In fact, we have two main well 
defined maxima at times t ~ 14.6 and t ~ 19.5. The peak in the 
reconnection speed is again caused by a fast outflow of gas visi- 
ble at time fi - 13.90 in Figure |8] (first column, top panel). If we 
look at the magnetic field topology at the same time (Figure [8] 
first column, middle panel) we see a magnetic island which ap- 
pears in the region of fast outflow. 

At next time step (f2 = 15.70) the topology of magnetic and 
velocity fields are still very mixed and complex (Figure [8] sec- 
ond column, top and middle panel), however the magnetic field 
does not exhibit the very strong bending seen in the previous 



?7„ = 5-10" ', Pi„i=0.5, k, .= 12 




8 20 22 



Fig. 7. The evolution of the reconnection speed for the model 
with Pi„j = 0.5, injection scale kf = 12 and uniform resistiv- 
ity rju - 5 ■ 10""*. Marked times correspond to those shown in 
Figure |8] 



step. As the simulation proceeds the reconnection rate slowly de- 
creases, reaching even slightly negative value^ (Figure [t] from 
f ~ 16tof ~ 18.5). 

At time step - 18.15 (Figure[8] last column) we show how 
our model looks when the measured reconnection rate is neg- 
ative. We can see a strong accumulation of the magnetic field 
around the midplane which extends over almost the whole com- 
putational domain. The region of velocity fluctuations is broad- 
ened in comparison to the previous time steps, even though the 
volume within which we drive turbulence is unchanged. This 
quite unstable situation sets off' the extreme growth of reconnec- 
tion rate (from t ~ 18.5 to f ~ 19.5) shown in Figure]?] The cur- 
rent density distribution is plotted in the bottom row in Figure [8] 



' This should not be taken too literally. Our measure of the recon- 
nection speed is misleading when the midplane magnetic field topology 
becomes progressively more tangled. 
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Fig. 6. Topology of velocity (top panels) and magnetic field (middle panels) in the presence of turbulence shown at three times 
fi = 12.15, t2 - 13.85, f3 - 15.15. In the lower panels we plot the absolute value of current density at those times. The variation 
of the reconnection speed Vj^ for this model is presented in Figure [s] Turbulence is injected with power Pi„j = 0.5 and at scale 
k - 12. The uniform resistivity rju is equal 5 ■ 10 The brighter shades correspond to the larger values of displayed quantities. 



K. Kulpa-Dybel et al.: Reconnecdon in weakly stochastic B-fields in 2D 




t,= 18.15 



1.2 1.0 



O.b 



0.6 >- 0.0 



-0.5 



-0.4 -0.2 0.0 0.2 0.4 |v^ 



t, = 13.90 



1.0 



0.5 



>- 0.0 



-0.5 



-0.4 -0.2 0.0 0.2 0.4 
X 



t, = 13.90 



0.5 




0.0 



170 



1 00 



80 



I 



-0.4 -0.2 0.0 0.2 0.4 




-0.5 



-0.4 -0.2 0.0 0.2 0.4 |v 

X 

t,= 1 5.70 



t,= 18.''5 



2.0 



1.5 0.5 




■.5 -0.: 



0.0 



2.0 1.0 



1.5 0.: 



.0 > 0.0 



0.5 



I 0.0 



2.0 



1 .5 



10.5 



lo.o 



.4 -0.2 0.0 0.2 0.4 




-0.4 -0.2 0.0 0.2 0.4 
X 

t,= 1 8.' 5 



-0.5 



-0.4 -0.2 0.0 0.2 0.4 




Fig. 8. Topology of velocity (top panels) and magnetic field (middle panels) in the presence of turbulence shown at three times 
fi = 13.90, t2 - 15.70, f3 = 18.15. In the lower panels we plot the absolute value of the current density for the same times. 
The evolution of the reconnection speed Vj^ for this model is shown in Figure [t] Turbulence is injected with Pi„j = 0.1 at the 
wavenumber k - 12. The uniform resistivity rju is equal 5 ■ 10 Brighter shades correspond to larger values of displayed quantities. 
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Fig. 9. Dependence of the reconnection rate on the power of in- 
jected turbulence P„y for models with and without turbulence, 
Vj^ and Vf, respectively. In all models we use the same values 
of the uniform resistivity 77,, = 5 ■ 10"^ and the injection scale 
ki„j - 12. The error bars indicate the variance of the reconnec- 
tion rate. In the case of Sweet-Parker reconnection (squares) the 
errors are neglected and not shown. 



Fig. 10. Dependence of the reconnection rate on the injection 
scale lj„j for models with and without turbulence, Vj^ and V^^, 
respectively. In all models we use the same values of the uniform 
resistivity 77,, = 5 ■ 10"^ and the power of injected turbulence 
f = 0.1. In the case of Sweet-Parker reconnection (squares) the 
variance is negligible and not shown. 



At every time step many small reconnection events occur, as in- 
dicated by local growth of the current density. 

3.2.1 . Dependence on the Strength of Turbulence 

In order to check how the reconnection rate Vj^ depends on the 
strength of turbulence we make several simulations with differ- 
ent values of turbulent power f ,„j (see model PD in Table[T]). The 
rest of input parameters have the same value in all these models. 

In Figure |9] we present the dependence of the averaged re- 
connection speed on the power of injected turbulence. Diamonds 
and squares correspond respectively to the reconnection rate 
with and without turbulence. Both of them are averaged over 
the fixed period of time: from t-llot-9 - laminar reconnec- 
tion (Sweet-Parker stage) and from t - 10 to f = 20- turbulent 
reconnection. The variance of the reconnection rate is calculated 
using the standard deviation method. In case of the Sweet-Parker 
reconnection the variance is negligible and is not shown here. 
On the other hand, in the presence of turbulence the reconnec- 
tion rate undergoes continuous strong variations and the variance 
is only slightly smaller than the mean reconnection rate. We find 
that the reconnection rate increases with power of turbulence and 
scales as ~ P^J^. . 

3.2.2. Dependence on the Injection Scale 

Similarly we determine the dependence of the reconnection rate 
on the wavenumber at which we inject turbulence kf. We 
compute several models with varying kj while keeping all other 
parameters the same (see model PD in Table [TJ. 

Figure 10 shows the resulting dependence. We can clearly 



see the strong relation between injection scale kf and the recon- 
nection speed. Namely, the reconnection rate and errors increase 
with decreasing injection wavenumber kf. Fitting to the simula- 



3.2.3. Dependence on the Resistivity 

The most important goal in our study is to check the dependence 
of the reconnection rate on the uniform resistivity. We run sev- 
eral simulations with different values of uniform resistivity rju 
and power of turbulence Pi„j (see model RD in Table [TJ. 

In Figure 1 1 we plot the obtained dependence of reconnec- 



tion rate on uniform resistivity. For Sweet-Parker configuration 
the averaged reconnection rate V^'' is represented by squares. 
Dispersion of calculated points is almost negligible and not 
shown. As we can see the reconnection speed V^^ increases with 
uniform resistivity. From the fitting we obtain that the reconnec- 
tion rate V^^ scales with the uniform resistivity as ~ tj^^. This 
is in agreement with theoretical prediction (Sweet 1958 , Parker] 
1957) and confirms that our 2D model works well in the laminar 
reconnection stage. 

Adding turbulence to the system leads to a weaker depen- 
dence between the reconnection rate and the uniform resistivity. 
Moreover, the power of turbulence also influences this relation- 
ship, as is clearly visible in Figure 



squares - P,-, 



11 



(triangles 



Pinj = 0.01 



0.1). We fit lines to calculated points and 



find that the dependence between reconnection rate and uni- 
form resistivity is stronger for lower values of P,„j. Namely, for 
Pi„j = 0.1 and Pi„j = 0.01 we get that the reconnection rate 

scales as ~ 77'/'' and ~ rjl/^, respectively. 

For low values of the uniform resistivity (rju < 1 ■ 10"^- 
Pi„i -0.1 and rju <1 ■ 10"'*-P,„/ = 0.01) we see that the obtained 
reconnection speeds Vj^ are almost the same. This similar val- 
ues of reconnection rate may be caused by prevailing numerical 
diffusion. Then, taking into account higher values of ?;„ does not 
modify the overall value of the total resistivity. However, in this 
case the numerical diffusion will also influence the reconnection 
rate in the Sweet-Parker configuration, what does not happen in 



our simulations (see Figure 1 1 



tion results gives Vf ~ 1^^ 



,2/3 

inj ■ 



Results described above indicate that in 2D turbulent recon- 
nection depends on the uniform resistivity. However, the ob- 
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Fig. 11. Dependence of the reconnection rate on the uni- 
form resistivity rju for models with and without turbulence, 
Vj^^ (diamonds and triangles) and V^^ (squares), respectively. 
Turbulence is injected at the scale k - \2 with P,„y = 0.1 and 



Fig. 13. Dependence of the reconnection rate on viscosity v for 
models with and without turbulence, Vj^ and ^, respectively. 
In all models we use the same values of the uniform resistivity 
riu - 5 ■ 10"^, turbulence strength P = 0.1 and the injection 



Pinj = 0.01 . In the case of Sweet-Parker reconnection (squares) wavenumber ki„j = 12. In the case of Sweet-Parker reconnection 



the variance is negligible and not shown. 



(squares) the variance is negligible and not shown. 
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Fig. 12. Dependence of the reconnection rate on the anomalous 
resistivity rja for models with and without turbulence, Vj^ (dia- 
monds) and V^'' (squares), respectively. Turbulence is injected at 
the wavenumber A-/ = 12 with = 0.1. The uniform resistivity 
J],, is equal 5 - 10 The dotted line represents the mean reconnec- 
tion rate obtained for the Sweet-Parker stage, dash-dotted marks 
the value of the reconnection rate without anomalous resistivity 
and dashed line corresponds to the mean reconnection rate cal- 
culated for four models with different values of the anomalous 
resistivity. 



tained dependencies are not certain because of the large disper- 
sion of results for small values of the uniform resistivity ?]„. 

We also test the dependence of the reconnection rate on 
the anomalous resistivity. We run several models with the same 
value of uniform resistivity ?/„ = 5 ■ lO"'* and critical current 
density jcm = 50 but for different values of the anomalous re- 



sistivity ria (see model AD in Table [TJ. In Figure 12 we plot 
the reconnection rate calculated for the Sweet-Parker configu- 
ration (squares) and in the presence of turbulence (diamonds). 
The dash-dot ted l ine determines the mean value of the recon- 
nection rate (Vj^f^^ ^q) - 0.051) obtained for models with the 
anomalous resistivity. As we see, the presence of anomalous re- 
sistivity causes an increase in the reconnection speed compared 
to a model with //« = - dotted line (Vj^ - 0.039). However, 
for different values of the anomalous resistivity the reconnection 
speed is almost the same. More precisely, we observe a nearly 
negligible increase in the reconnection rate and a clear increase 
in the variance at larger values of the anomalous resistivity rja. 

3.3. Dependence on the Viscosity 

Next we test the influence of viscosity on the reconnection rate 
for fixed values of the injection scale, turbulence strength and 
uniform resistivity (see model VD in Table [l]). The results are 
shown in Figure 13 where we plot averaged reconnection rates 
for the laminar (squares) and turbulent (diamonds) cases. In the 
Sweet-Parker configuration the reconnection rate V^'' is almost 
the same for viscosity v < 1 ■ lO"-' and starts to decrease for v > 
210"^. In the present of turbulence the fit to the calculated points 
shows that the reconnection rate Vj^ scales with the viscosity as 



-1/3 



Furthermore, when the viscosity is larger than 1 ■ 10 



-3 



the rate of reconnection for models with and without turbulence 
are similar. 



3.4. Tine Case of a Uniform Initiai Configuration 

In order to check the influence of our boundary conditions and 
method of driving turbulence on the reconnection process we 
compare two models with different initial magnetic field config- 
urations: antiparallel and uniform. The antiparallel initial con- 
figuration of the magnetic field is the same as in the case of the 
Sweet-Parker reconnection, the uniform one is horizontal with 
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Fig. 14. Comparison of the reconnection rates obtained for mod- 
els with uniform and antiparallel initial magnetic field configura- 
tions. All parameters describing models are the same: the power 
of turbulence P,„y = 0.1, the injection wavenumber k/ = 12, and 
the uniform resistivity coefficient r]„ = 5 ■ 10""*. The topology of 
the velocity and magnetic velocity fields, as well as the distribu- 
tion of the current density for both models are shown in Figure |6] 
(antiparallel initial magnetic field) and in Figure[T5](uniform ini- 
tial magnetic field). 



Bx - I- The rest of parameters characterizing both models are: 
77,, = 5 ■ 10""^, kf - 12, f = 0.1 and inflow/outflow boundary 
conditions. In both models we inject turbulence in the same way 
from f = 9 to f = 10. 



In Figure 14 we show the comparison of reconnection rates 
obtained for the uniform and antiparallel initial magnetic field 
configuration. In the first case the reconnection rate is almost 
equal zero. We can see some small fluctuations, however we 
cannot observe any significant growth of the reconnection rate. 
In the second case the value of reconnection rate is higher and 
changes violently. For more precisely description of this case see 
Sect. El] 

The almost negligible value of the reconnection rate for the 
model with a uniform initial magnetic field indicates that, as ex- 
pected, magnetic reconnection does not occur This also supports 
the conclusion that our calculation of the reconnection rate is not 
an artifact of our boundary conditions or our method for imple- 
menting turbulence. 

Figure [15] shows the topology of magnetic and velocity field 
as well as the distribution of current density for model with ini- 
tial uniform configuration of the magnetic field. As before tur- 
bulence is injected in a volume surrounding the midplane which 
causes perturbations of the velocity field in this region (see the 
left panel in Figure 



mean magnetic fie 



15 1. Most of fluctuations propagate along the 
d running more or less parallel to the mid- 
plane (see the middle panel in Figure 15 1. The well defined bend 
of magnetic lines presented in the case of turbulent reconnec- 
tion are not visible here. This is because we do not have a diffu- 
sion region, where the magnetic lines change their direction and 
strength, and can be bent by fluctuations. On the other hand, tur- 
bulence causes small perturbations of magnetic field which leads 
to an enhancement of the absolute value of the current density 
(see the middle panel in Figure 15 1. 



3.5. Long Term Variance of the Reconnection Rate 

The time evolution of the reconnection rate for the same model 
as in Figure |5]but calculated over a longer period of time (f(.„j - 
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Fig. 17. Time evolution of the total mass M (solid line), magnetic 
E„,„g (dashed line) and kinetic Eiti„ (dotted line) energies during 
Sweet-Parker (from f = to f = 9) and turbulent (from f = 9 to 
f = 50 reconnection with uniform resistivity 77,, = 5 ■ lO"'*. The 
kinetic energy Eki„ is not amplified in this plot. 



50) is shown in Figure 16 We see that even after 50 Alfenic 
times the reconnection rate does not reach a steady state. Strong 
and continuous fluctuations are seen over the whole simulation 
time. All these quasiperiodic changes of the reconnection rate 
are apparently driven by something like the tearing mode in- 
stability. Namely, big loops of magnetic field are continuously 
created (see Figure 18 1 and ejected through boundaries. 



In Figure 17 we show the evolution of mass, magnetic and 
kinetic energy during the Sweet-Parker (from f = to f = 9) 
and turbulent (from f = 9 to f = 50) reconnection. We see that 
adding turbulence to the system or using open boundary condi- 
tions does not introduce instabilities in the total mass or kinetic 
and magnetic energies. 

What is more, increasing tend from 20 to 50 only slightly 
influences the average reconnection rate and its variance. For 
models with f,„rf = 20 and te„d = 50 we get Vj'^ = 0.039 + 0.019 
and Vf^ = 0.041 ± 0.017, respectively. 

This indicates that extending the simulation time does not 
lead to a stable state of reconnection and does not change our 
results. 



4. Discussion 

4.1. The Goais of this Study 



The present paper provides a follow-up study to Kowal et al. 
( 2009[ l. In that paper fast magnetic reconnection of a 3D weakly 
stochastic magnetic field, predicted in the LV99 theoretical 
study, was confirmed. Can turbulent reconnection be fast in 2D 
as well? The question may sound rather scholastic, as the nature 
we deal with is definitely three dimensional. However, there are 
at least two reasons why answering this question is important. 
First of all, the claim in LV99 is that 3D effects are essential for 
fast reconnection. Thus it is important to test this claim and to 
explore whether fast reconnection in the presence of turbulence 
can be carried over to 2D. Then, in a number of earlier studies it 
was conjectured that magnetic reconnection could become fast 
in the presence of turbulence for purely 2D reconnection config- 
urations. 

In all respects apart from dimensionahty the present numeri- 
cal set-up is identical to that in |Kowal et al.| ( p009] l. In particular. 
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Fig. 15. Topology of velocity (left panel) and magnetic field (middle panel), and the absolute value of the current density distribution 
(right panel) in the case of the initially uniform magnetic field at the time t = 14. Brighter shades correspond to larger values of the 
displayed quantities. Turbulence is injected with the power P,„y = 0.1 at wavenumber kf = 12. The uniform resistivity j],, is equal 
5 ■ 10""*. For comparison, the model with the same set of parameters but with an initially antiparallel configuration of magnetic field 
is presented in Figure |6] 
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Fig. 16. Time evolution of the reconnection speed for the model with t P,„j - 0.1, an injection wavenumber kj = 12 and a uniform 
resistivity rj^ - 5 ■ lO""*. 



the excitation of turbulence in both case at subAlfvenic veloci- 
ties, thereby preventing field reversals arising from turbulence. 
For measurements of the reconnection rate we adopt the proce- 
dure presented in Kow al et al.| ( 2009^ . This enables us to make 
easy comparison between the 3D and 2D reconnection results. 



4.2. 2D versus 3D Reconnection 



Until very recently 2D geometry has been favored for magnetic 
reconnection studies. Its advantage stems from the fact that it al- 
lows one to achieve much higher Lundquist numbers compared 
to its 3D counterpart. The most important question is whether 
the nature of reconnection is the same in 3D and in 2D. It is 
suggestive that the answer to this question is positive if the re- 
connection of laminar magnetic field is involved. Our results on 



the Sweet-Parker reconnection above are very similaip]to the re- 
sults of the Sweet-Parker reconnection in |Kowal et al.| ( |2009l l. At 
the same time our comparison of the results in the present paper 
with those in Kow al et"aL] ( |2009] l shows that the reconnection in 
the presence of turbulence is very different. 

The aforementioned difference is not surprising if one takes 
into account that the nature of the 2D and 3D MHD turbulence is 
rather different. In particular, magnetic field wandering, which is 
an essential part of the LV99 reconnection model, has a radically 
different nature in 3D and 2D. What we have in 2D is coherent 
displacement of magnetic field lines by Alfvenic perturbations, 
while in 3D magnetic field lines can slip past one another, enjoy- 
ing the freedom provided by the additional dimension. As a re- 
sult, in LV99 different bundles of wandering magnetic field lines 

- The obtained dependence of reconnection pace on the uniform 
resistivity coincide in 2D and 3D and is consistent with theory, i.e. 
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Fig. 18. Topology of magnetic field in the presence of turbu- 
lence at time t - 49. The uniform resistivity is equal 5 ■ 10 
Turbulence is injected with the power P„y = 0.1 at scale k = 12. 
Brighter shades correspond to larger values of the displayed 
quantities. 



can enter the reconnection zone independently, and reconnect si- 
multaneously. It was shown in LV99 that it is this simultaneous 
reconnection of independent magnetic flux bundles that makes 
reconnection insensitive to resistivity. This effect is absent in 2D 
where the entry of fresh flux into the reconnection zone is con- 
strained by the rate of removal of reconnected flux. 

The second effect that enables fast reconnection in the LV99 
model is an increase of the thickness of the outflow region. 
There we see some similarity between 2D and 3D reconnection. 
Indeed, while the outflow in the Sweet-Parker reconnection is 
constrained by the thin slot determined by Ohmic diffusivity, it 
is due to an increase of the thickness of the outflow through tur- 
bulence that 2D turbulent reconnection is faster than the Sweet- 
Parker rate. The outflow region in both 2D and 3D case depends 
on the intensity of driving and the scale of turbulence injection. 
However, the functional dependencies are different and this dif- 
ference stems from the difference in the reconnection physics. 

As we mentioned earlier, magnetic field wandering plays 
a decisive role in 3D reconnection. In 2D reconnection is not 
fast, since we observe the creation of magnetic islands as a re- 
sult of turbulence. This process is rather limited unless the tur- 
bulent injection velocity approaches Va- We believe that this 
explains why the dependence of the reconnection rate on the 
Ohmic resistivity becomes steeper as the turbulence weakens 
(see Figure 111. Indeed, for weaker turbulence the alternate com- 
pression and expansion of the current sheet, which leads to mag- 
netic island formation is weaker and reconnection is similar to 
the Sweet-Parker model. Similarly, the combination of anoma- 
lous resistivity and magnetic island formation [preaches a satu- 
rated level of efficiency for modest values of the anomalous re- 



While the 3D reconnection of weakly stochastic magnetic field de- 
pends on the Alfvenic, i.e. the incompressible, component of magnetic 
turbulence, the creation of magnetic islands depends on the compress- 
ible component of the turbulence (see |Cho & Lazarian||2003] for the 
decomposition of MHD turbulent motions). 



sistivity, and appears to become insensitive to further increases 
(see Figure 12 1. 

It is important to note that magnetic field structures, i.e. mag- 
netic islands, that we observe in 2D simulations do not appear in 
3D sim ulation ( |Kowal et aL||2009i l. What is more, I Kowal et aT] 
( 2009 1 obtained a stable value of the reconnection speed in the 



presence of turbulence in a short period, about one Alven time. 
In our work we are unable to reach this state because the recon- 
nection rate continuously and violently fluctuates. Consequently 
the average speed of turbulent magnetic reconnection in 2D is 
burdened with much higher statistical errors than in 3D. 

In our work we determine the dependence of the reconnec- 
tion rate on four quantities: the power of turbulence P,„y, the in- 
jection wavenumber kf, the uniform ?/„ and the anomalous rja re- 



sistivities. The same analysis was done by Kowal et al. (2009 1 in 
three dimensions. For the power of turbulence and injection scale 
they found the following scalings: Vj^ ~ pI^J and Vj^ ~ i^j^j. 
Compared to our findings the dependencies obtained in 3D are 
stronger Namely, in the 2D case the reconnection rate grows 
with the power of turbulence as Vj 



TB 



tion scale as Vl^ ~ l 



a/3 
''"J ' 



1 /3 

p.. and with the iniec- 

mj J 



4.3. Earlier Studies of 2D Turbulent Reconnection and Their 
Relation to Work 

As we mentioned in the introduction, the concept of turbulent 
enhancement of reconnection is not unprecedented. In IKowal 



et al. (2009) we presented a list of papers where turbulent ef- 
fects were cited as the source of fast astrophysical reconnec- 
tion. Compared to the present paper, all these papers lacked a 
precise means of measuring the reconnection speed and there- 
fore the numerical simulations were providing mostly qualitative 
results. Moreover, they adopted periodic boundary conditions, 
which made it difficult to study turbulent reconnection for more 
than a single Alfven time. 

Among these papers, the pioneering work were the numeri- 
cal studies in Matthaeus & Lamkin|( |1985[|1986| ). There the anal- 
ysis of 2D simulations of turbulence revealed the formation of 
magnetic islands and X-points reminiscent of the Petscheck pro- 
cess. 

In our 2D model of turbulent reconnection we also observe 
continuous formation of magnetic islands which are ejected from 
the reconnection zone due to open boundary conditions. For in- 
stance, we observe that when big loop of magnetic field is re- 
moved from the box through boundaries, the reconnection rate 
increases sharply. After that the reconnection rate drops until 
new magnetic island can be created. 

The formation of islands is also a consequence of the tear- 
ing instability. In this vein, Loureiro et al. ( 2004 1 examined the 



nonlinear growth of tearing modes and obtained a fast growth of 
magnetic islands. However, in their simulations they used peri- 
odic boundary conditions, which keep such islands in the recon- 
nection zone. 



The most directly relevant work is Fan et al. (2004 2005 i. 
They found that the turbulent magnetic reconnection in solar at- 
mosphere could be described by a two phase process: first-slow 
and second-fast. The rapid stage of reconnection is caused by the 



coalescence instability ^Biskamp & Welter| ( |1980[ ); |Wu & Chang 



(2001 1 which may also enhance the reconnection rate in our sim- 
ulations. 

Here, and in contrast to the 3D model of LV99, we find that 
in 2D the influence of the uniform resistivity on the reconnec- 
tion rate is stronger for higher values of rju and lower values 
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of Pjnj. Similar results have also been obtained by (Loureiro et 



al.|2009 1. They claimed that the reconnection process is last m 



2D only in some particular circumstances i.e. for higher values 
of the Lundquist number and moderate levels of turbulence. We 
disagree with the claim of fast turbulent reconnection in 2D. In 
fact, in the most of models analyzed by ( Loureiro et al. 2009) the 
reconnection process depends on the uniform resistivity. Again 
an important limitation of the aforementioned study is their lim- 
ited averaging arising from their choice of boundary conditions. 

4.4. Implications 

Two dimensional magnetic reconnection is the result of a rather 
artificial configuration. The value of our present study is that 
it clarifies the role of the effects of dimensionality for the actual 
3D reconnection of astrophysical magnetic fields in the presence 
of weak turbulence. The latter has many essential implications 
starting from the First Order Fermi acceleration of energetic par- 



ticles ( de Gouveia dal Pino & Lazarian 2005 ; Lazarian & Opher 
|2009 ') to Solar Flares (see Lazarian et al. 2009 ) and removal of 
magnetic field during star formation (Lima et al. 2009) . These 
impUcations are discussed in more detail in Kowal et al. ( 2009| l. 



5. Conclusions 

In this paper we have examined the results of 2D simulations of 
reconnection process in presence of subAlfvenic turbulence. Our 
findings may be summarized as follows: 

- Unlike Sweet-Parker reconnection, the reconnection of a 
weakly stochastic magnetic field is fundamentally different 
in 2D and 3D, in agreement the LV99 study. Reconnection 
in 2D depends on resistivity and is not fast. 

- The enhancement of magnetic reconnection in 2D arises 
from an increase in the thickness of the plasma outflow layer 
due to the creation of magnetic islands which are ejected 
from the reconnection layer. 

- The power of turbulence P„y and the injection scale kf have 
an influence on the reconnection rate. In our study this is: 



V, 



TB 
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